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ABSTRACT 

0^ , We study the impact of correlated instrumental noise and non-circular antenna beam patterns on 

primordial non-Gaussianity analysis. The two systematic effects are reproduced in the case of the 
Planck mission, using PLANCK-like realistic simulations. The non-Gaussian analysis is conducted 
with different approaches. First we adopt a blind approach, using the Spherical Mexican Hat Wavelet 
. and the Minkowski functionals, and then a /nl estimator. We look respectively for false primordial 

non-Gaussian detections and for bias or variance increase in the estimated /nl value. Even if some 
slight effects are present, we can not observe any significant impact of the 1 // noise and the asymmetric 
beam on non-Gaussianity searches in the context of the Planck mission. 

Subject headings: cosmic microwave background — cosmology: observations — early universe — 
methods: data analysis — methods: statistical 
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1. INTRODUCTION 

Today the Gaussianity of the primordial fluctua- 
tions has become a key observable in cosmology. The 
standard single-field inflationary mod els predict unde- 
tecta b le deviation s from Gaussianity (|Acquaviva et al.1 
2003t iMaldacenal I2003D . while alternative scenarios 
can produce different levels of non-Gaussianity, usu- 
ally parametrized by the no n -linearity par a meter 
/W (IKomatsu fc Spergejl 120011: iBartolo et all 120011: 
Linde fc Mukhanovlll997tlLvth et al.ll2003l; IBabich et all 
20041: IChen et all I2007at IHqlman fc Tollevl 120081 
Chen et all l2007ri lLanglois et al.ll2008t) . Since different 
models predict different limits on the /nl parameter, 
the constrains on Gaussianity can discriminate between 
different scenarios. Nevertheless testing non-Gaussianity 
is challenging. Firstly because the primordial deviations 
from Gaussianity are expected to be tiny. Secondly, 
residual foregrounds and systematic effects could in- 
troduce non-Gaussian signatures which can mimic 
the primordial non-Gaussianity signal. Moreover it is 
interesting to look for any non-Gaussian signal, not 
only of the kind predicted by Inflationary models, 
but it does not exist a statistical test sensitive to all 
possible non-Gaussian manifestations of a random 
field. Therefore the problem must be approached with 
qualitatively different tests. The actual and forthcoming 
Cosmic Microwave Background (CMB) experiments 
have the capability to strongly constrain primordial non 
Gaussianity, thanks to their high sensitivity and angular 
resolution. In order to achieve this goal an accurate 
control on the systematic artifacts is essential. Today 
the best CMB data come from the Wilkinson Microwave 
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Anisotropy Probe {WMAP) satellite. The WMAP 
team analysis on the 5 year data found them consistent 
with Gaussianity, putting the following constrains on the 
non-linearity parameter / NL : -9 < f^l al < 111 (95% 

CL) a nd -151 < f^f 1 < 253 (95% CL) (IKomatsu et alJ 
2009). These parameters constrain non-Gaussian mod- 
els which predict large contributions to the bispectrum 
respectively fr om the squeezed and the equilateral 
config urations flBabich et al.l 1200^: iLinde fc Mukhanovj 
1997tlLvth et al.l 120031 : IBabich et all 120041 : IChen et all 



2007a; Hol man fc Tollevl 120081 : IChen et all 12007 



Langl ois et al.l 2008f) . These li mits have be e n imp roved 



with an optimal analysis by Smit h et"ail (|2009h and 
ISenatore et all (|20(59l) . finding -4 < fft£ al < 80 (95% 

CL) and -125 < f^ il < 435 (95% CL) respectively. 
Consistent results with different methods were fou nd by 
other groups fCurto et alJ 120091: iPietrobon eFal] 120091 : 
iRudiord et al.l 2009) . However other studies have shown 
some hints of non-Gaussian features at the 2 — 3er level. 
These include the presence of a very cold spot in the 
data, deviations of the expected statistics of peaks from 
Gaussian expectations, non-Gaussian distribution of the 
anisotropies in the southern emisphere a nd others (e.g. 
dVielva et all l2004t lHansen et alt l2009t iMcEwen etall 
2008)). Mor eover a detection o f a non -zero /nl has been 
achieved bv lYadav fc Wandeltl (|2008f ) on the WMAP 3 
year dataset. With the same statistical method adopted 
by the WMAP team, they found 27 < f^ al < 147 
(95% CL). The forthcoming data of European Space 
Agency (ESA) CMB mission Planck will help to clarify 
these puzzles. Planck will image the whole sky with 
outstandi ng sensitivity, angular resolu t ion an d frequency 
coverage (|The Planck Collaboration! 120051 ). To take 
advantage of this top-level performances, a fine control 
of all possible systematic effects and of their impact on 
the tests of non-Gaussianity is crucial. 

In this paper we consider two instrumental effects: the 
correlated noise, also known as 1// noise, and the asym- 
metry of the antenna beam. We have investigated their 
impact on primordial non-Gaussianity searches carried 
out with different statistical tests: Spherical Mexican 
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Hat Wavelets, Minkowski functionals an d a /nl estima- 
tor, o ften dubbed as "KSW" estimator (|Komatsu et alj 
2005). The systematic effects are studied in the context 
of the Planck mission. 

The paper is organized as follows. In Section [2 we 
describe the main characteristics of the Planck mission 
and the systematic effects under study. Then in Section[3] 
we shortly review the statistical tests used in this work. 
The analysis method applied and the results obtained 
are presented in Section 01 Conclusions are drawn in 
Section [5] 

2. THE SYSTEMATIC EFFECTS 
2.1. The Planck mission 

Planck, successfully launched on may 14, 2009, con- 
sists of two different instruments, called the Low and the 
High Frequency Instrument (LFI and HFI respectively). 
LFI covers the frequency range 30-70 GHz with an array 
of receivers based on High Electron Mobility Transistor 
(HEMT) amplifiers. The HFI receivers are instead based 
on bolometers, and operate between 100 and 857 GHz. 
Each feed horn collects radiation from the telescope and 
feeds it into one or more detectors. The satellite ob- 
serves the sky from the 2nd Earth-Sun Lagrangian point 
(L2), spinning at *~ 1 rpm with the spin axis kept in 
anti-Sun direction. The telescope field-of-view is offset 
from the spin axis by ~ 85°, so that the observed patch 
traces large circles on the sky. The detectors will make 
repeated observations of the same sky circles before re- 
pointing the spin axis. As the spin axis follows the Sun, 
the obser ved circles sweep through the sk y at a rate of 
Pday" 1 (|The Planck Collaborationl 120051 ). 

In order to achieve the precision aimed in the Planck 
CMB data, the increase in angular resolution and sensi- 
tivity must be supported by an excellent control of the 
systematic effects. In the following we describe two typi- 
cal instrumental effects: the 1// noise and the asymmet- 
ric beam. 

2.2. The 1// noise 

HEMT amplifier gain fluctuations, noise temperature 
fluctuations, thermal instabilities produce low-frequency 
noise in the data streams, also called 1/ / noise. In gen- 
eral the instrumental noise can be described as a sum 
of white and correlated noise, writing its Power Spectral 
Density in the form: 
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where fk is the "knee" frequency at which the white 
and the 1// noise contribution are equal, f s is the sam- 
pling frequency and a is the white noise standard de- 
viation per sample integration time (t — l// a ). The 
value of the spectral index a depends on the dominant 
physical process which generates low-frequency noise. If 
not corrected the 1/f noise leads, coupled with the ob- 
serving strategy, to unwanted stripe-wise structures in 
the final sky map. The effect of the 1/f noise can be 
app roximated by a con stant offset for each scan cir- 
cle (jJanssen et al.1 119961) ■ If the scanning strategy is 
redundant, as the Planck strategy actually is, it is 
possible to derive the amplitudes of these offsets con- 
sidering the intersections between different scan circles. 
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Fig. 1. — (top) Map of white and 1/f noise before destriping; 
(bottom) after destriping stripes are visible due to residual 1/f 
noise (the white noise contribution has been subtracted). We re- 
mark the different peak-to-peak values of the maps. 
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Fig. 2. — The black line is the power spectrum of a map with 
white and 1/f noise before (right) and after (left) destriping (for 
fk = 0.05 Hz, a = 1.7). The gray line is the power spectrum of a 
map with only white noise. Note the different values range of the 
two plots. 

This is the key fe ature of the d e stripi n g techniques for 
map-m a king (e.g. IMaino et al.l (|2002f ); [Ashdown et al.l 
(|2007al) ; IKeihanen et all (|2004D ) . 

A destriping method can signicantly reduce the 1/f 
noise. Anyway a residual contribution is still present 
after the correction. There are residual stripes in the 
map (Fig. [1]), and in the angular power spectrum we 
can observe an excess of noise respect to the white noise 
level at low multipoles I < 50 (Fig. [J). Analytical stud- 
ies have shown that the residual 1/f noise introduces a 
block-diagonal structure in the covariance matrix of the 
harmonic coefficients of the map. The off-diagonal terms 
have been found to be much smaller than the diagonal 
components, but not completely negligible in compari- 
son to the tiny n on-Gaussian signal we are looking for 
(Efstathiou 2005). Therefore is important to investigate 
the impact of these correlations on the non-Gaussianity 
statistical tests. 

2.3. The asymmetric beam 

Actual CMB observation requires multi-frequency fo- 
cal plane arrays to achieve proper angular resolution, 
sensitivity and spectral coverage. This is indeed the 
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Fig. 3. — Main beam contour plot of one of the Planck-LFI 
70 GHz feed horns (simulated assuming an ideal telescope). The 
dotted lines are the bivariate Gaussian fitting (color scale: -90-^0 
dB). 

case of Planck, where we have the HFI focal plane, 
with 36 horns, inserted into the ring formed by the 11 
LFI horns. As a consequence, some feed horn can not 
be located close to the focal plane center, so that the 
beam pattern is asymmetric rather than circular. Real- 
istic Planck optical simulations of the feed horns cou- 
pled to the telescope (jSandri et alj|2004t IBrossard et al.1 
2004) have shown that the main beam patterns are well 
fitted by an elliptical shape given by a bivariate Gaus- 
sian (Fig. [3]). Moreover the finite response in time of 
the detectors leads to an elongation of the beam along 
the scan direction. The horns scan the sky mainly in the 
same direction, perpendicular to the ecliptic plane. Only 
the ecliptic poles are swept in several directions, which 
makes the effective beam more symmetric. 

The convolution of the sky signal in a map T(f) with 
the beam can be written in harmonic space as 

m=e 

^) — ^ ^ ^ ^ Bi rn cic ni Y em (f), (2) 

i m=-l 

where Yt m are the spherical harmonics and Bg rn and ai m 
are the harmonic expansions of the beam and the sky 
signal respectively. For a circular symmetric beam we 
have Bi m = Bi rn i , but if the beam is asymmetric J5/ m 
is no longer isotropic and it depend s also on the az- 
imuthal direction m (|Wu et al.l l2001f ). The azimuthal 
dependence increases the deconvolution complexity and 
the computational cost, thus an exact deconvolution of 
an asymmetric beam is not straightforward. A beam 
not correctly measured and treated in the data analysis 
pipeline may artificially distort the underling structure 
of the CMB signal, introducing a spurious dependence 
on the azimuthal direction and messing up the isotropy 
of the field. We will investigate the im pact of this ef- 
fect in the non-Gaussianity searches (seelAshdow n et al.l 
(2009) for a study of the asymmetric beam impact on 
map-making) . 

3. STATISTICAL TESTS OF NON-GAUSSIANITY 

In this work we have made use of two qualitatively dif- 
ferent statistical approaches to test non-Gaussianity. In 



the first approach no assumption is made about the non- 
Gaussian model, and the Gaussian hypothesis is tested. 
In the following we will refer to such tests as blind tests. 
In this work we have chosen to use the Spherical Mexi- 
can Hat Wavelet (SMHW) and the Minkowsky function- 
al (MF hereafter). This tests make use of Monte Carlo 
(MC) simulations to calibrate the Gaussian behavior. A 
departure from the MC calibration is interpreted as a 
non-Gaussian detection. The blind approach has the ad- 
vantage to be model independent, but it has a weak sta- 
tistical power, i.e. it is difficult to interpret the mean- 
ing of an eventual detection of non-Gaussianity in the 
data. In order to put quantitative constraints on the non- 
Gaussianity level we have to adopt another approach. 
Testing specific non-Gaussianity arising from inflation- 
ary models we can put limits to the non-linearity pa- 
rameter /nl- We have used a cubic statistic to estimate 
the /nl parameter for both the local and the equilateral 
non-Gaussianity. 

Given the difference in the two approaches, also the 
possible impacts of the systematic effects arc different. 
For the blind tests the issue is in the MC calibration, 
since it is not computationally manageable to include 
the l/f noise or the asymmetric beam in sets usually 
of thousands of simulations. Therefore the presence of 
the systematics in the data can in general cause depar- 
tures from the MC calibration, that can lead to false 
non-Gaussian detections. Also in the non-blind cubic 
statistic it is often not feasible to take exactly into ac- 
count the systematic effects, due to the complexity of 
the numerical or analytical treatment involved. A not 
accurate handling of the systematics can potentially in- 
troduce a bias in the estimated /nl- Even if on average 
the effect is zero, and there is no bias, the spurious cor- 
relations introduced by the systematics can increase the 
variance of the estimator. An example of this behavior is 
the effect of the anisotropic noise, which introduces cor- 
relations between small and large scales. If not corrected, 
these leads to a noticeable increase of the variance at h igh 
multipoles (jCreminelli et alJl2006t iSmith et al.ll2009f ). 

Before analysing the impact of the systematics effects, 
we will briefly describe the statistical methods applied in 
the study. 

3.1. SMHW 

The main idea of the wavelet technique is to decom- 
pose the temperature field into a set of coefficients, which 
contain information about the signal at a given position 
and scale. We can perform a continuous analysis, tak- 
ing the scale and the position in a continuum of values. 
The wavelets are a convenient technique to study non- 
Gaussianity, since the wavelet transform is linear. Thus 
the wavelet coefficients preserve the Gaussianity prop- 
erties of the original field. Moreover they allow us to 
study the non-Gaussian features at different scales. In 
particular the SMHW enhances the non-Gaussian signa- 
tures with spherical symmetry. The SMHW is the exten- 
sion to the sphere of the Euclidian Mexican Hat W avelet 
(jCavon et al.ll2000HMartmez-Gonzalez et al.ll2002f l. It is 
given by 
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N(R)-: 



R 2 i? 4 

1 + T + T 



1/2 



where R is the scale, (9, <fi) represent polar coordinates 
on the sphere, (y = 2 tan | , </>) are polar coordinates in 
the tangent plane to the North pole, and N(R) is a nor- 
malisation constant. Given the temperature field on the 
sphere T(9, 0), the wavelet coefficients are defined as: 

w(B, <f>;R)= I d9' d(f>' sin 9' f (x + /i)*(6>'; R), (4) 

9' 

, sin </>), ii — 2 tan —(cos <p\ sin </>'). 

Several statistics of the wavelet coefficients Wi can be 
studied at a given scale R. In this work we have analyzed 
two simple non-Gaussian estimators: the skewness S(R), 
i.e. the degree of symmetry about the maximum, and 
the kurtosis K(R), which indicates the degree to which 
a statistical frequency curve is peaked. Assuming zero 
mean at each scale we have: 



Q 

x = 2 tan - (cos < 
2 V 



S(R): 



— Y 

Nr ^ 



1=1 
Nn 



Wj(R) 
a(R) 

Wj{R) 
a(R) 



(5) 



where cr 2 (R) — yJ2 i J t iWi(R) 2 J /Nr is the variance and 

Nr is the number of coefficients (equal to the number of 
pixels N p i X in our case). 

3.2. Minkowski junctionals 

A general theorem of integral geometry states that the 
morphology of a convex body in D-dimensional space can 
be completely des cribed by D + 1 Minkowski functionals 
(Minkowski 119031 ). In the case of the two-dimensional 
CMB field we consider the excursion sets Q, defined as 
the set consisting of the parts of the sky in which the 
temperature anisotropies amplitude u is above a thresh- 
old v. Q = Q{v) — {u\{u — (u))/ao > v}, where o 2 = 
(u 2 ) — (u) 2 . Each of the functio ns has an exp licit analytic 
expression for a Gaussian field flbmita 1986). Thus they 
provide an appropriate statistical test for Gaussianity. 
Different algorithms have been implemented to estimate 
the MF of a pixeled CMB sky map. In this work we 
have developed our algorithm following the approach of 
Eri ksen et all (|200l . using the derivative computation 
of the HEALPix 6 routines (|G6rski et al.ll2005f ). The first 
functional, called the area functional, is defined as: 



1 



dA 



N pix (u > v) 



bs JQ 



N 



tot 



(6) 



pix 



and is equal to the total area of Q. The second equality 
holds in the HEALPix pixelization, in which all pixels 
have exactly equal area. Then we define the length func- 
tional, i.e. the total length of the boundary dQ of the 
excursion set: 



C(u) = / dl 

JdQ 

6 http://healpix.jpl.nasa.gov 



(7) 



where C is normalized to have a unit integral value: 
/ L{y)dv = 1. Finally the third functional is the genus, 
which measures the connectivity of Q, being equal to the 
number of isolated regions (hot spots) minus the numbers 
of holes (cold spots). It can be defined by: 



S(u) = 



N max (u) + jV mil » ~ jVsad(^) 

A" max (-oo) + iV min (-oo) + N sad (-oo) ' 



(8) 



i.e. as the number of maxima and minima minus the 
number of saddle points in the excursion set, divided 
by the sum of all stationary points in the whole un- 
thresholded field. In addition we use another morpho- 
logical descriptor: the skeleton length, defined as the 
zero-contour line of the map S: 



S — T-^T-fiiT-^ — T-gg) + T-^g{T 2 e — T 2 ^), 



(9) 



where semicolons denote the covariant derivatives. We 
compute the length of the skeleton of that part of the 
field that lies above some threshold is, normalized by the 
length of the skeleton of the whole un-thresholded field. 
The zero-contour lines of S can be interpreted as the set 
of lines that extend from extremum to extremum along 
lines of maximum or minimum gradient. For details of 
the MF algorithm see lEriksen et alj (|2004f) . 

3.3. /nl estimator 

The three-point correlation, or the bispectrum, func- 
tion is the lowest-order statistic able to distinguish non- 
Gaussian from Gaussian perturbations. Indeed for Gaus- 
sian CMB the expectation value is exactly zero. There- 
fore we can look at the bispectrum as a natural observ- 
able to constrain the non-Gaussianity amplitude pre- 
dicted by inflationary models, determined by the /nl 
parameter. After measuring the bispectrum modes of a 
given map we can build a cubic estimator of /nl in the 
following way: 

/nl = (10) 
5 P rim is a cubic statistic defined as 



5, 



prim 



E 

ti<t 3 <i 3 



D 



obs 



p>prim 



Cg 2 Ci 3 



(11) 



where Bfff 2l is the bispectrum extracted from the ob- 
served map and Bj^ ai is the ansatz for /nl = 1 of the 
non-Gaussian model assumed, while N is a normalisation 
factor defined as 



N 



E 



tl<l2<l 3 



(12) 



where Ci is the theoretical power spectrum of the model, 
after including the instrumental beam and noise. 

We adopt the fast numerical imp lementation of 
the e sti mator in Eq. |JTU|) fou n d by [K omatsu et al.l 
<l2005h: iCreminelli et all (12006ft : lYadav et alJ ~if2008): 
iSmith & Zaldarriagal (|2006) . We consider two differ- 
ent primordial shapes of non-Gaussianity: local non 
Gaussianity, in which we expect the major contribu- 
tions to the bispectrum from squeezed configurations, 
where l\ <C ^2,^3; and equilateral non-Gaussianity, 
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which predict large bispectrum contributions for equilat- 
eral shapes, i.e. for l\ ~ I2 ~ £3. The corresponding /nl 
estimators are described in the literature cited above. In 
the following subsections we recall the definitions of the 
main quantities involved. 

3.3.1. local estimator 

In the context of a local model of primordial non- 
Gaussianity it is useful to define the quantities 

a,(r) = i J k 2 dkg Ti (k\ n (kr), (13) 

/3i(r) = ^ J k 2 dkP{k)g Ti {k)j e Xkr), (14) 

where gTi(f) is the radiation transfer function, ji(kr) is 
a Bessel function and P(k) is the power-spectrum of the 
primordial perturbations. With these quantities we can 
obtain two filtered maps from the harmonic coefficients 
at m of the observed CMB: 



At -\ oit,{r)Wi , 
A(r, n) = — -= atrnXtrnW , 



em 



C, 
Pe(r)W e 
C, 



(15) 
(16) 



where We is the beam function of the instrument and 
Ce is the theoretical power-spectrum Ce, after including 
effects of beam and noise as Ce = CeW 2 +a 2 . Here we are 
approximating the noise as homogeneous (n£ m n|, TO ,) ~ 
CnSee'Smm '- With these maps we can construct a cubic 
statistics ()Komatsu et al.l l2005) 



S'X = **j r 2 dr J ^A{r, h)B 2 (r, n). (1 



7) 



This expression reduces exactly to Eq. (fTTj) , so we can 
define the estimator of Eq. (|10p . In the case of full sky 
coverage and homogeneous noise, the estimator is fully 
optimal, i.e. it saturates the Cramer-Rao bound, provid- 
ing the minimum variance estimator. But it is not realis- 
tic to consider a full sky coverage. If f s k y is the fraction of 
sky observ ed we multiply the normalisation N of Eq. (fT2l) 
by fskyN (jKomatsu et all 120051 : ICreminelli et all [2006). 
Moreover in an experiment the noise results anisotropic, 
since different sky regions are observed for different 
amounts of time. Thus the noise covariance matrix 
(ne m rig, m ,} is no longer diagonal. The resulting corre- 
lation between small and large scales causes an increase 
of the variance of the estimator at high l's. This can be 
corrected with a linear t erm defined by (jCreminelli et ail 
12001 lYadav et~aT][200l 



^linear 



2 dr J d 2 h\2B(r, n){A slm (r, n) • 
(r,n)) MC + A(r,h)(B 2 lm (r,n)) 

MC }(18) 

where A S i m (r, n) and B S i m (r, n) are the A and B maps 
generated from Monte Carlo simulations that contain sig- 
nal and noise, and (. . .) denotes average over the simula- 
tions. The estimator is finally given by 

^local 1 ^linear 



This treatment greatly improves the final error bars of 
the estimator, but it do es not yet provide f ully optimal- 
ity This is achieved bv lSmith et al.l (|2009f ) through the 
implementation in Eq. (fTOf of the full covariance matrix. 



3.3.2. equilateral estimator 

In addition to ae(r) and (3e( r ) (Eqs. (I13I14|) ) we can 
define the functions 

lt {r) = - I k 2 dkP(k) 1/3 g Te (k)j e (kr), (20) 



6e(r) = -l k 2 dkP(k) 2/3 g T e(k)je(kr) . (21) 

and obtain the filtered maps C(r, n) and D(r, n) in the 
same way of the maps A(r, n) and B(r, n) (Eqs . (115I16|) ). 
Now the statistics is given by (|Creminelli e"t~al~l l2006) 



cquil 



prim 



-18 



r 2 dr 



d 2 h 



A(r, n)B 2 {r, n) 



2 o 

+ -£> 3 (r, n) - 2B(r, n)C(r, h)D(r, n) 
o 



(22) 



and the estimator results again f^f = S°^/N, where 
N is given by Eq. (|T2|) . In this case the anisotropic noise 
is not an issue, since the spurious effect is quite low for 
equilateral bispectrum configurations. 

4. ANALYSIS 

In order to test the impact of the 1/ / noise and the 
asymmetric beam we have analysed realistic Planck 
simulations of Gaussian CMB. The two systematic effects 
are studied separately, considering two different sets of 
PLANCK-like maps in which the 1// noise or the asym- 
metric beam are simulated. For each systematic the test 
consists in a comparison of the results obtained on the 
realistic simulation set, with those obtained on another 
set of simulations hereafter dubbed "nominal" , in which 
the considered systematic is not present. 

4.1. Simulation sets 

Realistic Planc k simulations can be obtained with the 
LevelS pipeline (jReinecke et al.ll2006D . an assembly of 
numerical tools specifically developed to model the out- 
put data of the Planck satellite. LevelS has been used 
to create simulations of the time-ordered-data (TOD) for 
the 12 radiometers of the LFI-70 GHz channel. The two 
systematic effects considered are expected to be stronger 
for the LFI channels. The choice in particular of the 70 
GHz channel in our study let us work with the best sen- 
sitivity and angular resolution among the LFI channels. 
The TODs cover 1 year of mission. Every set of TODs 
has been then converted i n a map, using the de stripcr 
map-making Springtide (jAshdown et al.ll2007aD . with 
N S ide = 1024 7 in the HEALPix pixelization scheme. For 
each systematic effect we have a set of 100 PLANCK-like 
maps. The size of the set is limited because of the large 
resources required by the TODs simulations and the sub- 
sequent map-making, both in terms of CPU computing 
time and disk space. We have to keep in mind this limited 
size when looking at the results. Nevertheless we have 
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found that 100 maps are sufficient to obtain meaningful 
results. We describe now the characteristics specific of 
each set. 

4.1.1. 1// noise simulation set 

For each radiometer the instrumental noise is simu- 
lated as a sum of white and l/f noise, with PSD given 
by Eq. ([T]). The parameters adopted are reported in 
Table [TJ Due to the Planck scanning strategy the 

TABLE 1 

PARAMETERS ADOPTED IN THE l/f NOISE SIMULATIONS 

standard deviation a = 1787.6 fiK 
knee frequency = 50 mHz 

l/f slope a = 1.7 

sampling frequency f s = 76.8 Hz 



noise is anisotropic. The white noise level chosen sat- 
isfies the goal sensitivity of the LFI-70 GHz channel 
(|The Planck Collaboration"! 12005( 1. Subsequent tests on 
the 70 GHz radiometers have shown a lower knee fre- 
quency. Moreover comparative map-making methods 
studies have found that the residual noise is slightly 
higher for Sprin gtide with resp ect to other destriping 
implementation (jAshdown et al.ll2007allbl . 120090 . There- 
fore these simulations are a conservative upper limit for 
the l/f noise. 

In addition to these 100 realistic noise maps, we have 
considered also a set of nominal simulations, i.e. 100 
white noise maps, in which the noise is a realization of 
anisotropic white noise with the same nominal standard 
deviation of the Planck simulations. To both noise map 
sets we have added 100 maps of Gaussian CMB sky. The 
CMB maps are random realizations of the WMAP-1 
year best fit power spectrum, convolved with the angu- 
lar resolution of the LFI-70 GHz channel with a circular 
Gaussian beam of 14'. Finally, for the SMHW and the 
MF analyses we need Monte Carlo (MC) simulations for 
calibration. Thus we have created 5000 MC maps of 
anisotropic white noise and Gaussian CMB. We remark 
that the MC calibration does not contain the l/f noise. 

4.1.2. asymmetric beam simulation set 

The LevelS pipeline has been used to simulate a 
PLANCK-like observation of 100 Gaussian CMB skies, 
created as random realizations of the WMAP-1 year 
best fit power spectrum. The sky observation is sim- 
ulated considering the best-fit elliptical beams of the 
LFI-7 GHz, obtained in realistic main beam simula- 
tions (jSandri et al.l 2004). The parameters of the bivari- 
ate Gaussian fitting for each radiometer are reported in 
Table [2] The twelve radiometers form symmetric couples 
due to the symmetry of the LFI focal plane. The sam- 
pling in time has been also simulated. Since the effect 
is a widening of the beam in the scan direction, it has 
been reproduced using a beam smeared in the rotation 
direction, perpendicular to the ecliptic plane. 

Furthermore we have created a set of 100 nominal sim- 
ulations of Gaussian CMB, made again as random real- 
izations of the WMAP-1 year best fit power spectrum, 
but in this case convolved with a circular Gaussian beam 
of 13'. This is the resolution of the symmetric beam 



TABLE 2 

PARAMETERS ADOPTED IN THE ASYMMETRIC BEAM SIMULATIONS 



id a 


FWHM" 


ellipticity c 


18a, 23a 


13.0328 


1.2601 


18b, 23b 


12.9916 


1.2184 


19a, 22a 


12.7053 


1.2509 


19b, 22b 


12.6494 


1.2162 


20a, 21a 


12.4795 


1.2460 


20b, 21b 


12.4273 


1.2224 



a LFI-70 GHz radiometers label; "a" and "b" refer to the polar- 
ization tuning. 

'Geometric mean of FWHM (full width at half maximum) 
of the major and minor axes of the beam ellipse (FWHM= 
V fwhm max fwhm min ) . 

c R,atio fwhm max /fwhm m ; n . 



6000 




Fig. 4. — Average power spectrum of 100 realistic simulations 
with asymmetric beam deconvolved with a circular Gaussian beam 
of 13' (black line). The theoretical CMB spectrum is shown with 
a dashed line. 

which best approximates the asymmetric one (Fig. |4|. 
For some analysis we have also added a random real- 
ization of anisotropic white noise to each CMB map of 
both sets. The nominal standard deviation per sample is 
er = 1787.6 /iK (as in Table Q}. We have finally created 
5000 MC simulations to calibrate the SMHW and the 
MF analyses, following the same recipe of the nominal 
simulations. 

4.2. SMHW analysis 

We have calculated the wavelets coefficients Wi at 18 
different scales R, starting from the pixel size: R = 3.4, 
6.9, 10, 13.7, 25, 50, 75, 100, 150, 200, 250, 300, 400, 
500, 600, 750, 900, and 1050 arcm. For each scale we 
have computed standard deviation a(R), skewness S(R) 
and kurtosis K(R) of the distribution of the coefficients 
(Eqs. (0)). For both the l/f noise and the asymmetric 
beam analysis, the statistics are calculated on the 100 
realistic simulations containing the systematic effect and 
on the 100 nominal simulations without it. For calibra- 
tion we have also calculated <r(R), S(R) and K(R) over 
5000 MC simulations. 

Since the distribution of the MC values is close to 
Gaussian, for every statistic and map we can construct 
a x 2 test in this form: 

X 2 (n) = (x(n) - (x MC )) T C^c (x(n) - (x MC » , (23) 

where the vector x has elements Xi = cr(Ri), S(Ri) or 
K(Ri) and C^c is the MC covariance matrix Cy = 
((xi — (xi)) (xj — (xj))). The x 2 test is performed on 
each map of the realistic and of the nominal simulation 
sets, with n = 1 — 100, and on each MC map, with 
n = 1 — 5000. Moreover we have established the accep- 
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tance intervals at different significance levels a for each 
statistic cr(R), S(R) or K(R) and scale. Fixing R, we 
look at the distribution of the values over the 5000 MC 
maps: an acceptance interval is defined as the one con- 
taining a percentage 1 — a of values, where the remaining 
percentage is the same above and below the interval, i.e. 
a/2 at each side. We have chosen three significance lev- 
els: 68% (la-), 5% (2a) and 1%. When a value falls out 
of an interval, the significance of the non-Gaussianity in- 
dication is given by a: lower is the significance level a, 
stronger is the indication. 

With these quantities we have constructed three figures 
of merit. Each figure is evaluated for standard deviation 
<j(R), skewness S(R) and kurtosis K(R) on both the re- 
alistic and the nominal simulation sets. We have looked 
at all the three statistics, but remembering that only 
skewness and kurtosis are Gaussian estimators. We have 
also studied the standard deviation for a better charac- 
terisation of the systematic effect. 

1. x 2 distributions comparison: 

We have tested if the 100 \ 2 values of the simu- 
lation set belong to the same distribution of 5000 
Xmc values °f the MC simulations. The degree of 
agreement is quantified as the percentage of Xmc 
greater than the mean of the x 2 of the 100 simula- 
tions; 

2. number of x 2 out of 2a: 

For each x 2 value of the simulation set we have 
looked at the percentage of MC simulations with a 
lower Xmc- ^ this percentage is less than 2.5% or 
more than 97.5%, then the % 2 value lies outside 2a 
of the Xmc distribution; 

3. number of outliers: 

For each scale we have counted the number of maps 
with statistic of the wavelet coefficients out of the 
acceptance intervals at 5% and 1%. 

The impact of the 1// noise or the asymmetric beam has 
been tested comparing the figures of merit of the realistic 
simulation set against the figures of the nominal one. 

4.2.1. SMHW: 1/f noise results 

The first figure of merit for the 1/f noise is showed in 
Fig. [51 The impact of the 1/f is noticeable only in the 
standard deviation: the x 2 distribution of the realistic 
set is not displayed in the plot because the values are 
out of range respect to the MC set (upper-left plot in 
Fig. [5]). This is indeed expected, because of the increase 
of the noise variance with respect to the white noise level 
of the MC calibration due to the 1/f noise contribution. 
All the other plots show an agreement between the set 
of maps and the MC calibration, indicated also by the 
percentage of Xmc g rea ter than the x 2 mean. Table [3] 
reports the second figure of merit. Apart again from the 
standard deviation, the 1/f noise does not increase the 
occurrence of skewness or kurtosis out of 2a from the 
Gaussian behavior as calibrated with MC simulations. 

The third figure of merit (Fig. [5]) shows the impact of 
the 1/f noise on the different scales. We can observe 
that at scales R ~ pixel scale the presence of 1// noise 
increase the probability for a map to have the standard 
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Fig. 5. — SMHW -1/f noise results: x 2 distributions comparison 
for realistic (top line) and nominal (bottom line) simulations. The 
dot-dashed line shows the (x 2 ) of the simulation set. We indicate 
the percentage of Xmc g rea t er than this mean. 



TABLE 3 

SMHW - 1/f NOISE RESULTS: NUMBER OF \ 2 0UT OF 2(7 



#of X 2 


std. deviation 


skewness 


kurtosis 


realistic 


100 


3 


3 


nominal 


1 


5 


7 



deviation of wavelet coefficients out of 1% or 5%. Ac- 
tually this scale is noise dominated and usually it is not 
considered in Gaussianity analysis. 

Before moving on in the analysis, a general comment on 
the third figure of merit is mandatory. We can sometime 
observe a number of outliers greater than what expected 
in particular for the nominal set, since these maps are 
obtained with the same pipeline of the MC simulations. 
This is only due to the limited size of the simulation sets, 
which increases the number of values in the distribution 
tails. 

We have repeated the analysis also with an amplified 
noise contribution: in both the realistic and the nominal 
simulation sets the noise maps have been multiplied by 
a factor equal to 2.5. The results found are qualitatively 
unchanged, thus we do not show them here. 

For a better understanding of the noise properties we 
finally have analysed only noise simulation, in which no 
CMB signal is added. We can look at the x 2 test of the 
realistic simulations in Fig. [7J Apart from the standard 
deviation, where again we have obtained values out of 
MC set range, only the kurtosis shows a somewhat lower 
percentage for the realistic set respect to the nominal 
one. Anyway the value still indicates a good agreement 
with the MC distribution. Also the number of x 2 out of 
2a reported in Tabled] is not affected by the 1/f noise 
contribution. But when we look at the third figure of 
merit in Fig.[8l we can see that the 1/f noise increases the 
number of outliers. The standard deviation is affected at 
all the scales, since here the increase of noise variance 
is no more hidden by any CMB signal. Also skewness 
and kurtosis are affected at some scales, but we observe 
that the number of outliers remains within statistically 
reasonable limits. 
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FlG. 6. — SMHW - 1/f noise results: number of outliers of the 
5% and 1% acceptance intervals. 
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FlG. 8. — SMHW - 1/f noise results: number of outliers [only 
noise simulations] 
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FlG. 7. — SMHW - 1/f noise results: x 2 distributions comparison 
[only noise simulations] for realistic (top line) and nominal (bot- 
tom line) simulations. The dot-dashed line shows the (x 2 ) of the 
simulation set. We indicate the percentage of Xmc greater than 
this mean. 



TABLE 4 

SMHW - 1/f NOISE RESULTS: NUMBER OF X 2 OUT OF 2cr [ONLY 

NOISE simulations! 



#of X 2 


std. deviation 


skewness 


kurtosis 


realistic 


100 


3 


4 


nominal 


4 


4 


7 



4.2.2. SMHW: asymmetric beam results 

The first figure of merit in Fig. [9] shows that the asym- 
metric beam has an impact on the standard deviation: 





39.32% 


1 
1 


i 

TJ 

L 



KURTOSIS 



10 20 . 30 40 50 60 



STD. DEVIATION 



..4.6.2.6%. . k 









1 




43.74% 


0- ] 






1 

n 


J 





KURTOSIS 
ftimlffi !ft - 



X 



10 20 . 30 40 50 60 



FlG. 9. — SMHW - asymmetric beam results: x 2 distributions 
comparison for realistic (top line) and nominal (bottom line) sim- 
ulations. The dot-dashed line shows the (x 2 ) of the simulation set. 
We indicate the percentage of Xmc g rea ter than this mean. 

the x 2 values of the realistic set are out of range with 
respect to the MC set values. We can see in Fig. @]that 
indeed the elliptical beam used in the realistic simulation 
produces a different smoothing with respect to the circu- 
lar beam, used in the nominal and in the MC simulations. 
Instead skewness and kurtosis percentages indicate a still 
good agreement with the MC calibration, even if lower 
with respect to the symmetric beam maps. The number 
of x 2 out of 2(7 in Table [5] confirms that the asymmet- 
ric beam has a significant impact on the x 2 distribution 
of the standard deviation statistic only. Looking at 
the number of outliers in Fig. [10] we observe that the 
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TABLE 5 

SMHW - ASYMMETRIC BEAM RESULTS: NUMBER OF \ 2 °UT OF 2a 





std. deviation 


skewness 


kurtosis 


realistic 


100 


3 
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nominal 
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Fig. 10. — SMHW - asymmetric beam results: number of outliers 

asymmetric beam increases the occurrence of wavelet co- 
efficients statistics out of the acceptance intervals. This 
happens only at small scales for the standard deviation, 
while we have a general slight increase for kurtosis and 
skewness. Anyway for these two statistics the numbers 
of outliers are still admissible and the departure from the 
MC Gaussian calibration is not significant. We have re- 
peated the analysis adding an anisotropic white noise to 
the CMB maps, both in the realistic and in the nomi- 
nal simulation sets. We have found that the noise does 
not modify the impact of the asymmetric beam and the 
results do not change. 



Gaussianity of the values distribution over the MC, we 
have constructed a % 2 test. We have performed a di- 
agonal x 2 test, due to the narrow step-size used in the 
thresholds range, defining 



x 2 H = E 



f{v,n) - (/mcM) 



^McM 



(24) 



where / is one of the MF, (fMc{ v )) and <Jmc{ v ) are the 
mean and the standard deviation over the MC simula- 
tions. The x 2 test has been performed on the realistic 
and the nominal simulation sets and on the MC calibra- 
tion maps. 

The evaluation of the results is achieved analysing two 
figures of merit, identical to the first two figures consid- 
ered in the SMHW analysis in i)4.21 i.e.: 

1. x 2 distributions comparison; 

2. number of x 2 out of 2a. 

The figures of merit are evaluated for each smoothing 
scale and functional. For both the 1// noise and the 
asymmetric beam study, we are going to compare the 
results of the figures obtained on the realistic and on the 
nominal simulation set, i.e. respectively with or without 
the systematic considered. 

4.3.1. MF: 1/f noise results 
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Fig. 11. — MF -1/f noise results: \ 2 distributions comparison 



4.3. MF analysis 

In order to investigate different angular scales, for both 
the 1// noise and the asymmetric beam study we have 
smoothed the sets of 100 simulations, realistic and nomi- 
nal, with Gaussian beams at {i.e. no smoothing), 10, 30, 
60, 120, 180 and 240 arcm. For every map and smoothing 
scale we then have estimated the MF, i.e. area, length, 
genus and skeleton length. The derivative computation 
is performed up to an appropriate multipolc l max for 
each scale: from 850 to 60 for the larger beam. The MF 
are estimated for 200 threshold values v = — 4<j -j- 4<7, 
but in the subsequent analysis we consider only thresh- 
olds —3(7 < v < 3a (—2.5(7 < v < 2.5a for the two 
larger scales). We have applied the same procedure to 
the 5000 MC maps for calibration. Having verified the 
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Fig. 12. — MF - 1/f noise results: number of \ 2 out of 2a. 
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In the first figure of merit (Fig. [TTj) it is straightfor- 
ward to see that there is not any noticeable difference 
between the results obtained with the 1// noise included 
and the results from the nominal simulations, in none of 
the smoothing scales applied to the maps. Also looking 
at the number of x 2 out of 2a in Fig. we can say that 
the agreement of the simulation set with the MC calibra- 
tion is not affected by the 1// noise. We have checked 
the robustness of the results repeating the analysis on 
simulations with an higher noise level. The noise maps, 
both in the nominal and in the realistic simulation set, 
have been multiplied by a factor of 2.5. The results found 
are almost identical to that showed. 

4.3.2. MF: asymmetric beam results 

Comparing the distribution of \ 2 values of the simu- 
lation sets with the MC calibrations (Fig. [T3| . we ob- 
serve that the asymmetric beam of the realistic simu- 
lations does not cause a lower agreement with respect 
to the nominal simulations. In Fig. Q3]we can see that 
for the asymmetric beam maps in many cases we have 
found more % 2 values out of 2a respect to the symmetric 
beam maps, but the increase is not remarkable. We have 
also analysed the same CMB simulations after adding 
anisotropic white noise, but the results do not change 
significantly. 
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MF - asymmetric beam results: \ distributions com- 
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Fig. 14. — MF - asymmetric beam results: number of x 2 out of 
2<r 
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Fig. 16.— /S uil - 1// noise results 



4.4. /ml estimator analysis 

In order to test the impact of the 1// noise and the 
asymmetric beam on the /nl parameter, we have per- 
formed the estimation with the cubic statistic described 
in £ 13.31 on the realistic and the nominal simulation sets, 
in both the local and the equilateral configuration. In the 
local configuration, when the noise is added, we have also 
corrected for correlations introduced by the anisotropic 
noise, computing the linear term of Eq. fll8[) . The bis- 
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Fig. 18. — ff m ' - asymmetric beam results 



pectrum is calculated up to different maximum values of 
multipole i ma x ■ We have proceeded looking at the mean 
of the /nl values obtained, to check against biases with 
respect to the expected null value. Then we have investi- 
gated the impact of the systematic effects on the variance 
of the estimator, comparing the standard deviations of 
the two sets of 100 /nl values obtained from the realistic 
and the nominal simulations. 



As explained in Q'3.31 correlation properties of the noise 
random field have to be properly accounted for, introduc- 
ing correction terms, like the linear term of Eq. (|18|) ; or 
computing the full noise covariance matrix and perform- 
ing an inverse cov ariance matrix weighting of the data 
( Smi th et alll2009D . But in the case of the 1// noise the 
estimation of the noise covariance matrix is numerically 
challenging, in par ticular at the high reso lution reached 
with Planck fsee IKeskitalo et alJ ([2009) for the study 
of the noise covariance matrix for Planck low-resolution 
data analysis). Before making any attempt in this direc- 
tion, it is therefore useful to evaluate the effective impact 
of the 1// noise without accounting for it in the /nl es- 
timation. Actually the noise is considered in the compu- 
tation of the Fisher matrix N (Eq. [H?)) and of the maps 
A(r,h), B(r,n), C(r,n) andD(r,n) (Eqs. {151 flS))). but 
it is accounted by means of its power spectrum, there- 
fore without considering the correlations introduced by 
the 1// noise in the realistic simulations. 

In Fig. [TS] we show the results for the local configura- 
tion, while in Fig. \W\ there are the equilateral configura- 
tion results. In both the cases we can not observe any 
bias in the mean of /nl- This is indeed expected, since 
the 1// noise is correlated but still Gaussian. But it is 
important to check against possible biases introduced by 
the not proper noise treatment itself. Moreover the 1// 
noise does not increase the /nl standard deviation, and 
the error bars have the same amplitude obtained from the 
nominal simulations. The small deviations between the 
obtained values and the limit given by the Fisher matrix 
are compatible due to the limited size of the simulation 
sets. 

4.4.2. /nl-' asymmetric beam results 
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We can repeat here the same considerations made for 
the 1// noise. The asymmetric beam can introduce cor- 
relations in the data, but the numerical description of 
this effect is demanding. Therefore we want to evalu- 
ate the effective impact of this systematic when it is not 
properly accounted for. In the estimation of /nl the 
beam considered in the computation is circular. In par- 
ticular we have considered a circular Gaussian beam of 
13', the same used in the nominal simulation. 

The results obtained for the local and the equilateral 
/nl are shown in Fig. [TTland Fig. IT51 respectively. Look- 
ing at the mean of the /nl values, not any bias is observ- 
able with respect to the null value in both the configu- 
rations. For the equilateral case, the standard deviation 
is also not affected by the asymmetric beam. In contrast 
in the local configuration we can observe an increase of 
the standard deviation at high multipoles. 

In order to understand if this increase is due only to the 
different smoothing of the elliptical beam with respect to 
the circular one, we have applied a transfer function in 
the /nl estimation. This function is calculated as the 
mean of the spectra of the realistic simulation set, di- 
vided by the theoretical spectrum convolved with the 
circular 13' beam. Therefore the "effective" beam, ob- 
tained as the circular beam multiplied for this transfer 
function, has the same amplitude of the elliptical beam. 
We have used this "effective" beam in the /nl estima- 
tion of the asymmetric beam maps. The results obtained 
using the transfer function are shown in green color in 
the same Figs. [17] and [18] We can observe that, even if 
lower, the increase in the standard deviation of the local 
configuration is still present. 

Then we have analyzed the same simulations after 
adding an anisotropic white noise contribution to the 
CMB maps. This is indeed a more realistic situation. 
The noise does not change the results for the equilateral 
configuration. On the other hand the standard devia- 
tions of the local /nl are significantly improved with re- 
spect to the situation without noise. We show that the 
local configuration results in Fig. [19] the increase of the 
error bars at high multipoles is no longer observable. 

5. CONCLUSIONS 

We have studied the impact of the 1// noise and the 
asymmetric beam on searches for non-Gaussian primor- 
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dial signals. These two systematics effects are indeed ex- 
pected to be present in the forthcoming Planck CMB 
data. With realistic PLANCK-like simulations we have 
evaluated if they can affect the non-Gaussianity analy- 
sis carried out with different statistical approach: blind 
tests, as the SMHW and the Minkowski functionals, and 
an /nl estimator for both local and equilateral configu- 
rations, known as the "KSW" estimator. The two sys- 
tematics have been studied separately. In the case of the 
blind tests we have looked at different figures of merit, 
checking for false non-Gaussian detections, while for the 
/nl estimator we have checked for the bias in the es- 
timated value and for the increase in the variance. In 
each case the results on realistic simulations have been 
compared with the results on simulations without the 
systematic effect examined. 

For both the 1// noise and the asymmetric beam anal- 
ysis, no tests conducted shows a significant impact on 
the final results. The 1// noise can slightly affect only 
the SMHW analysis, and only when the noise dominates 
over the CMB signal, e.g. at pixel scale. The blind tests 
are not significantly affected by the non circular beam, 
but in the case of the /nl estimator we can observe an 
increase of the variance of the local /nl at high multi- 
poles, when we consider simulations without noise. Any- 
way this effect is canceled once we consider more realistic 
noisy simulations. 

In conclusion, assuming that the simulations analysed 
well reproduce the actual Planck observations, if any 
non-Gaussian signal will be detected in the forthcoming 
Planck CMB data with these statistical tests, it will not 
be due to the 1/f noise or the asymmetry of the beam. 
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